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ABSTRACT 

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). 
Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for 
random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 
50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster 
size and porosity using these aggregate models. We find that while the shape of the phase function depends 
mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter 
and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to 
previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic 
aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2irci(>/\ ~ 1.6 
where ao is the monomer radius. Previous claims for such monomer size effects are in fact the combined 
effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase 
function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including 
the negative polarization observed for comets at scattering angles 160° < 9 < 180°. These aggregates have 
moderate porosities, V w 0.6, and are of sub-/im-size for the debris disk case, or /im-size for the comet case. 



Subject headings: dust, extinction - polarization - scattering - circumstellar matter - 
medium - stars: individual (AU Mic, GJ 803) 



comets - interplanetary 



1. INTRODUCTION 

Interplanetary dust particles (IDPs) c ollected in the 
Earth's strato s phere by high-flying aircraft (lBrownleelll985l : 
IWarren et al.l |T994) usually have irregular shapes and 
fluffy structures. Similar structures have been pro- 
duced in laboratory and microgravity experiments of dust 
partic l e interactions dWurm & Bluml 119981: iBlum & Wurml 
I2000t iKrause & Bluml |2004l) . it has also been sug- 
gested that interstellar dust grain s may consist primarily 
of such aggregate structure s (e.g.. iMathis & Whiffed [19891 : 
iDorschner & Henningl 19951) . with a mixture of various chem- 
ical compositions and vacuum. 

Porous, composite aggregates are often modeled as a clus- 
ter of small spheres ("spherules" or "monomers"), assembled 
under various aggregation rules. The optical properties of 
these aggregates can be calculated using numerical schemes 
such as the generalized mu ltisphere Mie (GMM) solution 
dMackowskil 1 1 99 ll : IXul 1 1 9971) or the discrete dipole approx- 
imation (PDA) metho d (e.g., iPurcell & Pen nyp acker! 1 19731 : 
iDraine & Flataul["l994f) . These methods have been used to 
study the optical properties o f different kinds of aggregates 
during the past decade (e . g.. IWesj 11991b ILumme & R ahola 
1994; Petrova et all 120001: iKimura et al]l2006t iBertini et alJ 
2007t lLasue et alJ l2009): most of those studies are dedicated 
to interpret the phase function and polarization of light scat- 
tered by cometary dust. 

In a companion paper (IShen et alJl2008l hereafter Paper I), 
we constructed aggregates using three specific aggregation 
rules: ballistic agglomeration (BA), ballistic agglomeration 
with one migration (BAM1) and ballistic agglomeration with 
two migrations (BAM2). We developed a set of parameters to 
characterize the irregular structure of these aggregates. While 
the BA clusters are essentially the Ballistic Particle-Cluster 
Agglomeration (BPCA) clusters frequently used in the liter- 
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2009), the newly-introduced BAM1 and BAM2 clusters have 
geometries that are random but substantially less "fluffy" than 
the BA clusters. The effective porosity V (eq. 12 in Paper I) 
increases from BAM2^BAM1^BA and covers a wide dy- 
namical range, allowing us to investigate the effects of poros- 
ity on the optical properties of the aggregates in a systematic 
way. Using these aggregation rules, we can construct grain 
models with various sizes and compositions. In Paper I, we 
computed total scattering and absorption cross sections for the 
three types of aggregates (BA, BAM1 and BAM2), for three 
different compositions (50% silicate and 50% graphite; 50% 
silica te and 50% amorphous carbon AC 1 . iRouleau & M artin 
1199 ll and 100% silicate), and for wavelengths from 0.1 /im to 
4 /jm. The purpose of this paper is to investigate the detailed 
light scattering properties of these aggregates, i.e., the phase 
function and the linear polarization. 

The paper is organized as follows: in $2] we recapitulate 
our aggregate models; the scattering phase function and lin- 
ear polarization for various ballistic aggregates are presented 
in §|3] where we explore the dependence of light scattering 
properties on aggregate properties; we present examples of 
aggregates that can be applied to circumstellar debris disks 
and cometary dust in §@] and we show that moderate poros- 
ity aggregates can reproduce the observed scattering and po- 
larization properties of dust in both solar system comets and 
extrasolar debris disks. We summarize our results in $5] 

2. AGGREGATE MODELS 

A detailed description of the target generation algorithms 
and resulting geometric properties of the BA, BAM1 and 
BAM2 clusters can be found in Paper I. Here we review some 
of the basic concepts that will be used in the following sec- 
tions. 
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TABLE 1 
Cluster Geometries 



cluster 


V 






"dip 


Figs. 


BA.256.1 


0.8598 


1.9249 


115656 


451.8 


6 


BA.256.2 


0.8525 


1.8925 


105855 


413.5 


6 


BA.256.3 


0.8553 


1.9050 


102725 


401.3 


6 


BAM1.32.4 


0.6188 


1.3791 


13903 


434.5 


4ab 


BAM1.256.1 


0.7060 


1.5038 


103921 


AOS Q 
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BAM1. 256.2 


0.7412 


1.5693 


107160 


418.6 


6 


BAM1.256.3 


0.6980 


1.4904 


107047 


418.2 


6 


BAM2.256.1 


0.5632 


1.3179 


113696 


444.1 


l,3,4cd,5,6,7 


BAM2.256.2 


0.5781 


1.3333 


103509 


404.3 


l,3,4cd,5,6,7 


BAM2.256.3 


0.5818 


1.3372 


107332 


419.3 


l,3,4cd,5,6,7 


BAM2.512.14 


0.6127 


1.3719 


211211 


412.5 


4ab 


BAM2.1024.1 


0.6386 


1.4039 


414977 


405.3 


4cd,8 


BAM2. 1024.2 


0.6476 


1.4158 


412077 


402.4 


4cd,8 


BAM2. 1024.3 


0.6387 


1.4041 


430184 


420.1 


4cd,8 


NOTE. — Naming convention: "BA.256.1" 
the N = 256 BA clusters. 
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Each aggregate is composed of N spherical monomers with 
radius ao- We define the "effective radius" of a cluster, a e ff, 
to be the radius of an equal-volume solid sphere; thus our 
aggregates have 



a e ff = N 1/3 a 



(1) 



The structure of the cluster is characterized by a porosity 
parameter V (see eq. 12 of Paper I) and a characteristic radius 
R = a e ff/(l — V) x ^ (see eq. 11 of Paper I), which depends on 
V and is typically 1-2 times a^tf. Tables 1 and 2 of Paper I 
give tabulated mean values of V and /?/a e ff for the three types 
of aggregates with 2 3 <N<2 lb . For a given value of N, the 
BA clusters have the highest V, while the BAM2 clusters have 
the lowest V . Information for the specific cluster geometries 
employed in this paper can be found in Table Q] including the 
porosity V, the number A/di P of dipoles representing the real- 
ization, and the number of dipoles per sphere, «dip = N& p /N. 
The actual geometry (including images) of these and other re- 
alizations of BA, BAM1, and BAM2 clusters can be obtained 
online. 1 

In Paper I we considered three different compositions: 50% 
silicate + 50% graphite, 50% silicate + 50% AC1, and 100% 
silicate. Silicate material accounts for perhaps 2/3 of the total 
mass of interstellar dust, and it is natural to assume that sil- 
icates will also provide the bulk of the refractory material in 
comets or debris disks. Interstellar silicates are amorphous; 
and amorphous silicates are believed to dominate the silicate 
mass even in the case of comets or circumstellar disks where 
crystalline silicates have been detected. We use the "astrosi li- 
cate" dielectric function dDraine & Ledll984lbraindl2003l) . 

Carbonaceous material provides a significant fraction of the 
total mass of interstellar grains, and this may also be true 
of dust in comets and debris disks. The smallest carbona- 
ceous particles in the ISM consist primarily of polycyclic 
aromatic hydrocarbon material, but the form of the carbon 
in the larger gra ins (where most of the carbon re sides) re- 
mains uncertain: Pendlet on & Allamandolal d2002l) conclude 
that the hydrocarbon material is ^85% aromatic (rin g-like) 
and 15% aliphatic (chain-like), but Dartois et al. (2004) claim 
that aliphatic material predominates, with at most 15% of the 
carbon in aromatic form. To explore the effect of material 
that is strongly absorptive in the visible, we use the dielectric 
tensor of graphite for the carbon in our mixed-composition 
aggregates. 

1 http://www.astro.princeton.edu/~draine/agglom.html 



In the present paper we study the scattering properties of 
aggregates with two compositions: 100% silicate, or 50% sil- 
icate + 50% graphite (volume fractions). In Paper I we found 
that aggregates consisting of 50% silicate + 50% AC1 amor- 
phous carbon had scattering properties intermediate between 
the 100% silicate and 50% silicate + 50% graphite aggregates. 

Calculations are per formed using DDSCAT version 7.0 
dDraine & Flataul 120081) . DD SCAT is a code based on the 
discrete dipole approx imation dPurcell & Pen nvpackerl l 19731; 
iDraine & Flataul Tl994T) . designed to compute scattering and 
absorption of electromagnetic waves by targets with arbitrary 
geometry and composition, for targets that are not too large 
compared to the wavelength A. For each cluster type (defined 
by N, aggregation rule, and composition) we generally av- 
erage over three random realizations and 54 orientations for 
each realization. 2 DDSCAT 7.0 allows us to treat the graphite 
monomers as randomly-oriented spheres with the anisotropic 
dielectric tensor of graphite. 

3. SCATTERING PHASE FUNCTION AND POLARIZATION 

The phase function and linear polarization of the scat- 
tered light as functions of scattering angle 9 can be re- 
trieve d from the elements of t he 4x4 Muller matrix Sy 
(e.g., iBohren & Huffman! 119831) . For unpolarized incident 
light, the scattered light phase function is proportional to 
Su = (47r 2 / A 2 )afC sca /aff2 (where dC^/dVl is the differential 
scattering cross section for unpolarized incident light) and 
the linear polarization parameter is p = — S^i/Sn. By defini- 
tion, the polarization is perpendicular/parallel to the scattering 
plane when p is positive/negative. 

We have obtained S\\(fi) and p{9) for our realization- 
and orientation-averaged aggregates for wavelengths 0.1 < 
A/jtfln < 4. For illustrative purposes, in most cases we will 
present the results for the BAM2 aggregates - the aggregate 
geometry with the lowest porosity. Orientation-averaged scat- 
tering properties for the clusters studied in this paper (includ- 
ing wavelengths not shown in the figures) are available on- 
line. 3 

3.1. Wavelength Dependence 

We first show the wavelength dependence of Su and lin- 
ear polarization for N = 256-monomer BAM2 clusters with 
monomer radius ao = 0.02 /im in Figure Q] for selected wave- 
lengths. For the N = 256 BAM2 case, a e ff = 0.127 pm, the 
porosity V « 0.58 (see Table [TJ, and the characteristic radius 
R w 1 .334a e ff = 0.17 fim . The phase function shows a rela- 
tively smooth dependence on wavelength A: for A < 0.5 /im 
(x = 2-kR/X > 2), it shows a strong peak in the forward scat- 
tering and a mild backscattering enhancement, with the over- 
all forward-backward asymmetry decreasing monotonically 
as the incident wavelength increases. For linear polariza- 
tion, the situation is more complicated. The polarization near 
9 « 90° first decreases as A increases, reaches a minimum at 
A w R, and then rises again with increasing A, approaching 
p(90°) = 100% in the Rayleigh limit. The wavelength A m i n po i 
where p(90°) is minimum is well-defined for the pure sili- 
cate, with Amin.poi w 0.17/im. For the graphite-silicate com- 
position it is less well-defined, with minima near ~ 0.17/xm 

2 9 values of the angle between the cluster principal axis Si and x (the 
direction of the incident light), and 6 values of the rotation angle /3 of the 
cluster around aj . We use a single value of the rotation angle <J> of ai around 
x because we average over 4 scattering planes. 

3 http://www.astro.princeton.edu/~draine/SDJ09.html 
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and ~ 0.45 /mi. The increase of polarization with increas- 
ing wavelength in the optical band is known as the polariza- 
tion color effect in cometary scattered light observations (e.g . , 
IChernova et al.|[l996t iLevasseur-Re gourd & Had ancikll2001l) . 
The reverse behavior of increasing polarization with decreas- 
ing wavelength in the UV band, however, is more complicated 
to interpret. It could be caused by the change in the size pa- 
rameter x = 2ttR/X, or changes in the dielectric function as A 
varies, or both. We will return to this point in 33.31 

It was shown in Paper I that the EMT-Mie model provides 
a good approximation for the total extinction cross section 
as a function of A, provided that the vacuum fraction / vac is 
set to / vac « V. We now test to see if the EMT-Mie model 
reproduces the scattering phase function Su(9) and polariza- 
tion p(9). For the EMT-Mie calculations we use an optimal 
value of vacuum fraction / vac = 0.55 and the same amount of 
solid material as in the N = 256 BAM2 clusters. For the ef- 
fect ive dielectric permittivity e e ff we use the Bruggeman rule 
(see lBohren & Huffmanll 19831) 

^ e, + 2e eff 
i 

where /; and e,- are the volume fraction and dielectric permit- 
tivity of each composition, including vacuum. There is al- 
ways only one solution of e e g that is physically meaningful. 
For graphite, we make the usual \— \ approximation, and take 

e = e(E || c) for / = ±/ grap hite, e = e(E _L c) for / = |/ gra phite- 

The EMT-Mie results are shown in Fig. [2] for the silicate- 
graphite and the pure silicate cases, in parallel to Fig. Q] 
The EMT calculations at fixed wavelength show resonances 
that arise from the use of spheres, but these should be 
smoothed out when modeling nonspherical particles, which 
are randomly-oriented and will not show such well-defined 
resonances. Therefore we have smoothed the EMT results us- 
ing a Gaussian kernel 



J aTnaexp 




ln(a /a)] 2 /2a 2 


Sij(a) 


J oTnaexp 


-[ln(a/a)] 2 /2 


a 2 





(3) 



where a = a/(a+ X/2ir) and a = (1 — /vac) _1 ' 3 fleff is the radius 
of the Mie sphere. The phase function and polarization are 
then computed using the smoothed Sy. 

By directly comparing Fig.[Tjand Fig. |2]it is evident that the 
EMT-Mie results for Sn(9) and linear polarization p(9) do 
share the same trends we see in the DDA calculations. Nev- 
ertheless, there are substantial differences between the EMT- 
Mie results and our DDA results. One obvious feature is that 
the EMT-Mie model tends to underestimate the backscatter- 
ing for short wavelength (A < 0.6 /im, x > 1.8), a feature al- 
ready revealed by the behavior of the asymmetry parameter 
g = (co&9) discussed in Paper I (fig. 12). For example, con- 
sider the forward-backward asymmetry 5'ii(0)/5'ii(180 o ) for 
A = 0.168 /im: the DDA calculations for the mixed graphite- 
silicate BAM2 cluster give ~ 450, while the EMT-Mie calcu- 
lation gives ~ 7000. 

To compare the EMT-Mie results and our DDA calcula- 
tions in detail we plot the relative differences in Fig. [3] for 
the silicate plus graphite case (upper) and the pure silicate 
case (bottom). The difference can be substantial for specific 
wavelengths or scattering angles. For example, for 9 w 90° 
scattering at A « 0.50/im, the EMT-Mie calculation underes- 
timates the polarization by a factor ~ 2, for both the graphite- 
silicate clusters and the pure silicate clusters. Although Paper 



I showed that EMT-Mie calculations can be used to obtain 
moderately accurate total extinction and scattering cross sec- 
tions, Figure 3 shows that the scattering phase function and 
polarization estimated using EMT-Mie calculations do not ac- 
curately reproduce the scattering properties of irregular clus- 
ters. 

3.2. Does Monomer Size Matter? 

There is another parameter that might affect Su an d po- 
larization: the monomer size. There have been claims that 
large monomer size is crucial in decreasing the polarization 
and in producing the negative polarization branch observed i n 
cometary dust (e.g., iPetrova et aT1l2000t iBertini et ai]|2007l) . 
However, in these previous studies, variations of monomer 
size were always coupled with changes in porosity V and 
cluster size R, hence effects attributed to varying the monomer 
size may in fact be due to variations in V or R. We have al- 
ready seen in Paper I that the apparent effects of monomer size 
on total cross sections are essentially the effects of varying V 
orR. 

To isolate the effect of monomer size, we compare clusters 
with the same a e ff and very similar V (thus R is also com- 
parable), but different monomer size ao- Thus the effect of 
monomer size, if there is any, is decoupled from other effects. 
We first consider the same example used in figure 8 of Paper I: 
the N =32 BAM1 cluster realization BAM1.32.4 (V = 0.619, 
R/ciefi = 1-379) with monomer size ao = 0.0504 /im and the 
N = 512 BAM2 cluster realization BAM2.512.14 (V = 0.613, 
R/ciefi = 1.372) with ao = 0.02 /im. Both clusters have a e ff = 
0.160 /mi and R = 0.220 /im. The orientation-averaged re- 
sults are shown in Fig. |4^,b for two wavelengths and for the 
silicate-graphite composition only. Although there are slight 
differences, the two cases have similar phase functions and 
polarizations: at constant R and V, varying the monomer size 
ao had little effect on the phase function and polarization. 

The above example employed moderate-sized clusters (x = 
2-kR/X < 3.9) composed of small monomers (27rao/A < 0.9). 
Fig- Et'd compares the scattering properties of 2 large clus- 
ters (R « 1.4/im, x = 11.1 and 13.9) with similar porosities 
V k, 0.6 but different monomer sizes. For A = 0.631 /im and 
0.794 /tm, clusters with ao = 0.10/mi and 0.16 /mi show sim- 
ilar (though not identical because of t he slight difference in 
porosity of the two clusters; see {sl3.41 i polarization p(9), de- 
spite the substantial difference in monomer size. 

3.3. Dependence on Cluster Size 

As we have discussed in ^3. II for fixed-size clusters, the 
dependence of the phase function and linear polarization on 
wavelength A is likely caused by the changes in both the 
size parameter and dielectric function. To investigate the ef- 
fects of cluster size at fixed incident wavelength (i.e., the ef- 
fects of size parameter alone), we use N = 256, BAM2 clus- 
ters with monomer size ao = 0.02, 0.025, 0.03 /im, or R ~ 
0.169, 0.212, 0.254 /im. These clusters have the same poros- 
ity, and as argued in the previous section, monomer size has 
negligible effects, hence any difference must be caused by 
changes in R. The results are shown in Fig. |5]for both compo- 
sitions. We present results at two wavelengths: A = 0.126 /im 
(< R) and A = 0.631 /mi (> R). In both cases the back- 
ward/forward scattering asymmetry increases with increasing 
the size parameter 2wR/X. 

In general, we expect p(90°) — > 1 in the Rayleigh scatter- 
ing limit R/ A <C 1 , with the peak polarization decreasing with 
increasing R/X. This decline with increasing R/X is seen in 
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Figure ]5b,d. However, the results in Figure [5k,c show that 
the variation of p mdx with increasing R is not monotonic: at 
A = 0.126 /mi when R > A and when the dielectric function 
is very absorptive, for both the 100% silicate and 50% sili- 
cate + 50% graphite N = 256 BAM2 clusters, the polarization 
is an increasing function of R over the range 1.3 < R/X < 2. 
Thus for these cases the polarization at A = 0.126/im has a 
minimum at some size /? cr ; t < 1.3 A. 

However, the dependence of polarization on R/ A depends 
on the dielectric function (and therefore on both composition 
and wavelength). For A = 0.631 /im the 100% silicate BAM2 
clusters with V « 0.6 have the polarization declining with in- 
creasing R out to R/ A = 2.2 (the largest value computed, see 
34.21 and Figs. [TJl and [8) - without showing a reversal in the 
polarization behavior. The situation is even more complicated 
for the silicate plus graphite case, where there is no coherent 
trend when R ks A (see Figs. [TJi, [7^ and|7J)). Based on the 
cases investigated thus far, it appears that when the dielectric 
function has only weak absorption (e.g., 100% astrosilicate 
at A = 0.631 /im), for fixed porosity V the polarization is a 
monotonically decreasing function of cluster size R from the 
Rayleigh limit 7?<cAupto,R/A<2. On the other hand, when 
the dielectric function is strongly absorptive (e.g., materials 
at A = 0.126 /im or silicate-graphite clusters at A = 0.631 /im), 
for fixed porosity the polarization declines with increasing R 
from the Rayleigh limit until it reaches a local minimum at 
R ss A (the transition is less distinct for the silicate-graphite 
case than for the pure silicate case), and then rises as R is 
further increased, at least out toS/A~2 (e.g., Figs. [SJi and 
EJi,b). 

3.4. Dependence on Porosity 

We now investigate the effect of porosity on the scattering 
phase function and linear polarization. Previous studies on 
the porosity effect using only the BPCA and/or the even more 
porous "ballistic cluster-cluster agglomeration (BCCA)" clus- 
ters were quite limited in the dynamical range of porosity, 
and changes in porosity were coupled with changes in clus- 
ter size. To decouple from the cluster size effect we choose 
clusters with comparable sizes, but different porosities. We 
use N = 256 BA, BAM l o and BAM2 clusters, with monomer 
size ao = 143, 170, 200 A respectively; hence these clusters 
have comparable size R ~ 0.17/im, but different porosities 
V = 0.85, 0.74, 0.58. 

We consider two regimes: A < R and A > R. The results are 
shown in Fig. |6]for two example wavelengths, A = 0.126 /im 
and A = 0.447 /im, and for the two compositions. It is evident 
that in both regimes, porosity has little effect on the shape of 
the phase function 4 . On the other hand, higher porosity tends 
to increase the linear polarization for both R > A and R < A. 

Most of the cases shown in this section have large 
linear polarization fraction [p(90°) > 50%]. Cometary 
dust typically has p(90°) < 40% and a negative branch 
of polarization at scattering angle ~ 160° - 180° (e.g., 
iLevass eur-Reg ourd & Hadancikl [20 01), obse rved at optical 
wavelengths. From the results of 33.U33.4I we expect that, 
in general, a reduced peak polarization and appearance of a 
negative polarization branch for 160 < 9 < 180° can be ob- 
tained by (1) increasing the cluster size R and (2) making the 
cluster more compact (lower porosity V). Examples of ballis- 
tic aggregates that are able to reproduce these cometary dust 

4 The slight difference in 5n/Sn(0) is likely caused by the different ge- 
ometry of the BA, BAM1 and BAM2 clusters. 



features will be presented in g4T 

4. APPLICATIONS OF BALLISTIC AGGREGATES 

The light scattering properties of our ballistic aggregates 
can be applied to various observations. Here we focus 
on debris disks and cometary dust, where single scattering 
dominates in the optically thin regime. We show exam- 
ples of aggregates that can reproduce qualitative and quan- 
titative features observed in the scattered light fr om the de- 
bris disk around A U Mic dGraham et al.l |2007|) and from 
cometary dust (e.g., iLumme & Rahola 1994; Petrovaetal. 



I2000t IKimura et al. 2006t lLasue et al.ll2009l) . Due to com- 



putational limits, we cannot probe a sufficiently large param- 
eter space to claim that our models are unique; nor do we 
attempt to fit a sophisticated model to the observations of a 
specific comet. Nevertheless, our examples (in particular the 
moderate-porosity BAM2 clusters) nicely reproduce most of 
the features observed in light scattered by debris-disk dust and 
cometary dust. 

4.1. Debris disk around AU Mic 

Polarization maps of the debris disk surround ing the nearby 
M star AU Microscopii have been obtained by iGraham et al.1 
(2007) using HST ACS in the F606W optical band (A c = 
0.590 /im, AA = 0.230 /im). The sca ttered light is polar- 
ized perpendicular to the disk plane. Graham et al. (2007) 
adopted the form for the phase function introduced by 
iHenvev & Greensteinl d 19401) : 



1 



1- 



11 47r(l+g 2 -2gcos(9) 3 /2 ' 
assumed that the polarization vs. scattering angle varies as 



(4) 



on 



sin 



1+ cos 2 9 



(5) 



and simultaneously fitted the phase function and linear po- 
larization as function of scattering angle to the observational 
data, obtaining g S3 0.68 and p mdx w 0.53. Graha met alj 
d2007l) suggest that very porous {V « 0.91 -0.94) /im-sized 
spherical grains or aggregates can produce these features 
based on Mie theor y and DDA calculations for BA clusters 
(IKimura et al.ll2006l) . 

Here we will show that random aggregates with a much 
lower porosity (V S3 0.6) can, in fact, better fit the observa- 
tions of the AU Mic debris disk. To reproduce the observed 
features, we require that the phase function and linear polar- 
ization are both close to the functions (0]i and © with the th e 
best-fit values of g and p mRX found by Graham etail (|2007). 
In particular, the intensity of scattered light at 9 = 0° should 
be approximately a factor of 150 larger than the intensity at 
9 = 180°, and the maximum polarization should be R3 0.5, al- 
though p(9) need not necessarily peak at 9 = 90°. 

Dust grains in debris disks will have a distribution of sizes. 
Much of the interstellar grain mass can be approximated by 
a power-law size distribution dn / dR cx R~ a for R < 0.25 /im 
with a R3 3.5 (Mathi set al.ll 1977b . Size distributions with a » 
3.5 can be obtained fro m models with c oagulation and col- 
lisional fragmentation (lDohnanyilll969t iTanaka et alJ Il996fc 
Weidenschillin^ll997t). 

For modeling comets and debris disks, we will consider a 
size distribution dn/dR oc /?" 3 5 . For a fixed porosity (i.e., a 
particular type of aggregate with a fixed N), this size distribu- 
tion is just dn/dao oc d^ 5 , where ao is the monomer size. 
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Hence the averaged phase function Sniff) and polarization 
are: 



Sniff) 



J ™* daoidn /dao)S\ \ (ao , ff) 



P(ff)-- 



J" m " dao(dn /dao) 
f"™" daoidn / dao)S\ i (ao , 0)p(ao, ff) 
da (dn /dao)S u (a , ff) 



(6) 
(7) 



where a m ; n and a max are the minimum and maximum values 
of monomer size in our size distribution. 

We consider N = 256 BAM2 clusters (V = 0.58), with a = 
0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08 /im, which correspond 
to R w 0.169, 0.254, 0.339, 0.424, 0.508, 0.593, 0.678/im. 
We show the calculated phase function and polarization for 
each of these clusters (averaged over orientations and 3 re- 
alizations) in Fig. [7] at two wavelengths A = 0.501 and 
0.631 jitm. For comparison with the scattering properties in- 
ferred for the dust around AU Mic, we calculate Sniff) an d po- 
larization piff) averaged over a size distribution dn jdR oc R~ 35 
with fl m ; n = 0.015 pm and a max = 0.065 /im (j? m j n = 0.127 /mi, 
^max = 0.551 /im), which are shown as dashed lines. The 
best-fit Henyey-Greenstein p hase function © and polariza- 
tion fitting function Q from Grah am et alJ ( |2007l) are shown 
as solid black lines. We plot the comparison for both the 
silicate-graphite composition (upper panels) and the pure sil- 
icate composition (bottom panels). 

As we can see from Fig. [JJ the size-averaged silicate- 
graphite clusters produce close matches to the Henyey- 
Greenstein model at these optical wavelengths for both the 
phase function and linear polarization. These clusters have 
porosity V w 0.6 and overall size R ~ 0.2- 0.5 /im, i.e., they 
are sub-/im-sized clusters with moderate porosity. 

As discussed in 33.41 increasing porosity will increase the 
polarization of scattered light. Our highest-porosity clusters 
are those BA clusters, which are commonly used in the liter- 
ature, referred to as "ballisti c part i cle-cluster agglomeration 
(BPCA)" c l usters 



OssenkopH 



Lasue et al. 



1993 



( e.g., [W est 199l llKozasa etal.ll 19921 



1993 



2007 



iKimura et all 120061: iBertini et all 
2009b . We found that if we replace the BAM2 
clusters in Fig. [JJ with BA clusters with the same number 
of monomers and monomer sizes, we can reproduce simi- 
lar phase function featur es but over-p r edict t he polarization. 
This is already seen in Grah am et al.l d2007) (e.g., their fig. 
8). Thus we conclude that compact BAM2 clusters fit the ob- 
servations of AU Mic better than the more porous BA clusters 
that have been considered previously. 

4.2. Cometary Dust 

The phase function and linear polarization of scattered light 
have been observed in a variety of comets. Though differ- 
ent comets show quantitative differences in the phase function 
and polarization, there are some common features: 

• The phase function shows strong forward scattering 
with a weak enhancement in the backscattering; the ge- 
ometric albedo [defined as A = (5 n [180°]A ? )/(47rG) 
where G « nR 2 is the averaged geometric cross sec- 
tion o f the grain] of backscatter ed light is less than 0.06 
(e.g. jHanner & Newburnlll989l) . 

• The linear polarization p(0) is a bell-shaped curve 
as function of scattering a ngle, with typical max- 
imum value of 10-30% (Dob rovolskv et all 119861: 



ILevasseur-Regourd et al.|[T996l) . although gas contam- 
ination in polarimetric measurements with wide-band 
filters might depol arize the observed scattered light 
dKiselev et alj|2004l) . 

• Within the 4000-7000A window, the polarization 
increases with wavelength, w hich is the so - called 
polarization color effect (e.g., jChernova et all 119961: 
ILevasseur-Regourd & Ha dancik 200ll). 

• Many comets show a negative branch of polarization at 
scattering angle larger than 150- 160°, with a minimum 
of > -2% (e.g. JPollfus et al.lll988l;rEaton et al.|[l992b . 

Most of these features, in particular the negative polariza- 
tion branch, have been successfully reproduced using vari- 
ous aggregates which differ in geome t ry, composition and 
porosity (e.g . . iLumme & Raholal Il994t 
IKimura et all l2006t IBertini et al.l l2007t 



'etrova et a 



Lasue et al 



2000; 
2009). 



The aggregates studied here are capable of producing all these 
features as well. In addition, we have demonstrated the ef- 
fects of grain size and porosity, and pointed out that the 
monomer size effect claimed by previous authors is in fact 
due to cha nges in cluster size R an d/or porosity V. Fo r ex- 
ample, in IPetrova et all ([2000) and IBertini et all (120071) . the 
difference of the prominence of the negative branch is caused 
by the effect of grain size when they increase the monomer 
size for the same configuration/porosity. 

The reason that those authors did not find a negative branch 
of polarization for small monomer size (ao < 0. 1 /im) and 
moderate number of monomers (a few tens) is that computa- 
tional limits prevented them from using a sufficient number 
of monomers. To test this, we have computed a few realiza- 
tions of BAM2 clusters with N = 1024, composed of small 
monomers («o = 0.08 /im). The results are shown in Fig. 
[8] for optical wavelength A = 0.631 /im, where the negative 
branch is evident for both compositions (it is more promi- 
nent at the usually observed red band A = 0.55 /im). In Fig. 
[8] we have also computed for the N = 1024 BAM2 clusters 
using monomer size ao = 0.1 /im, shown as open squares; 
the dashed lines are the average of the results of the two 
monomer sizes, for a size distribution dn/dR oc R~ i 5 running 
from R m i„ = 0.98 /im to R max = 1 .54 /im. 

These clusters are /im-size grains, whic h is consistent with 
the values found by other authors (e.g.. ILumme & Raholal 
11994 IPetrova et al.ll2000b IKimura et al Jl2006l) . although the 
exact values depend on composition as well as porosity. 
We find that the N = 1024, a = 0.08,0.1 /im BAM2 clus- 
ters composed of silicate and graphite (Fig. [8^) with R « 
1 . 1 - 1 .4 /im and P ss 0.6, reproduce a backscattering albedo 
A = [Sn(180 o )A 2 ]/(47T 2 fl 2 ) ~ 0.04 and a small negative polar- 
ization for scattering angles ~ 155 - 180° peaking at (~-l%), 
representative of the typical values found in cometary ob- 
servations, although the maximum polarization (~ 45%) is 
a little higher than observed. We may need somewhat more 
compact aggregates, or different composition (e.g., more sil- 
icate, see the right panel of Fig. [8) to lower the peak polar- 
ization. Alternatively, one may consider t he mixture of fluff y 
aggregates and compact solid grains (e.g.. lLasueeTaIl l2009). 
In their study, a larger fraction of porous BCCA aggregates 
is needed to produce a higher peak polarization for comet 
C/1995 Ol Hale-Bopp than for comet lP/Halley, consistent 
with our argument that high porosity helps increase the polar- 
ization. Since only very porous BCCA clusters were used 
in their modeling, it will be interesting to see if our more 
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compact BAM1 and BAM2 clusters will provide better fits 
for these comets in modeling the mixture of fluffy aggregates 
and compact solid grains. 

5. CONCLUSIONS 

We have studied the phase function and linear polariza- 
tion properties of light scattering by ballistic aggregates. We 
studied the wavelength dependence, cluster size dependence, 
and porosity dependence of the light scattering properties us- 
ing the discrete-dipole-approximation, and compared with the 
EMT-Mie model. Our main conclusions are: 

1 . It is shown that though the EMT-Mie model reproduces 
similar trends in these dependences, it differs quanti- 
tatively from the DDA calculations. We recommend 
using DDA calculations if accurate results are desired. 

2. Monomer size has negligible effects on the scattered 
light properties as long as monomers are small com- 
pared with the incident wavelength A. Even when the 
monomers are no longer small (e.g., 2nao/X ~ 1.6), 
monomer size appears to be of secondary importance 
for the phase function and polarization p{9). 

3. The phase function is mainly determined by R/X; in- 
creasing R/X decreases the backscattering relative to 
forward scattering. 

4. When R/X <C 1 (e.g., in the Rayleigh limit), increasing 
R/X decreases the polarization. For R > A, the depen- 
dence of polarization on R/X depends on the dielec- 
tric function: for materials that are not strongly absorb- 
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ing, increasing R/X results in decreasing polarization, 
at least for R/X < 2 (e.g., Figure |7J;,d, showing 100% 
silicate BAM2 clusters at A = 0.501 /im and 0.631 /im); 
however, at vacuum UV wavelengths where the materi- 
als are strongly absorbing, increasing R/X can increase 
the polarization (e.g., Figure|5Ji,c, showing scattering at 
A = 0.126 /mi). 

5. The degree of polarization depends on the size parame- 
ter as well as porosity, but high porosity helps increase 
polarization in both the R < A and > A regimes. 

6. We present aggregates with BAM2 geometry, moderate 
porosity V « 0.6, and sub-/im sizes which can repro- 
duce the scattered light phase function and polarization 
observed in the AU Mic debris disk. 

7. We present aggregate models with BAM2 geometry 
and moderate porosity V w 0.6 that can reproduce the 
albedo and polarization p(9) observed for cometary 
dust, including the negative polarization observed at 
scattering angles 160° < 9 < 180°. These aggregates 
are composed of silicate and graphite, and are of > /im 
size. Such moderately porous aggregates are promising 
candidates for cometary dust. 
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FIG . 1 . — Wavelength dependence of Su [normalized using Si i (0)] and polarization for the N = 256 and uq = 0.02 /im BAM2 clusters for two compositions, 
averaged over 3 realizations (BAM2.256. 1-3) and 54 random orientations for each realization. Upper: the silicate-graphite case; bottom: the 100% silicate case. 
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FIG. 2. — Same as Fig.[T] but for the EMT-Mie results (smoothed as in eq.[5J computed for / vac = 0.55 and same amount of solid material as in the N = 256 
BAM2 clusters in Fig.Q] 
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FIG. 4. — Tests for the effect of monomer size, (a) Comparable phase function and polarization for two silicate-graphite clusters - BAM 1.32. 4 and BAM2.512.14 
-with similar porosities V ~ 0.619, 0.613 and cluster size R ~ 0.221. 0.220 /im, but different monomer size ciq = 0.0504. 0.02 fim, for A = 0.355 fim (x = 2itR/X = 
3.9). (b) Same as (a), but for A = 0.631 fim (x = 2.2). (c) Comparable phase function and polarization at A = 0.631 fim for two large silicate clusters with similar 
porosity V ~ 0.57 and 0.64, and size R ~ 1.4/zm (x = 13.9), but different monomer size ag = 0.10,0.16/mi (2ttoo/A = 1.0, 1.6). (d) Same as (c), but for 
A = 0.794 fim. The clusters in (c) and (d) are each represented by 3 realizations (BAM2.256.1-3 and BAM2. 1024. 1-3) and 54 orientations per realization. These 
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FIG. 5. — Cluster size dependence of Sn and polarization for N = 256 BAM2 clusters with ag = 0.02,0.025,0.03 (im, for two compositions, averaged over 3 
realizations (BAM2.256. 1-3) and 54 random orientations. These clusters have typical sizes R ~ 0.169, 0.212, 0.254 /im, but same porosity V ~ 0.58. Examples 
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FIG. 6. — Effect of porosity on Si i and polarization for N = 256 BA, BAM1 and BAM2 clusters which have the same cluster size R ~ 0. 17 fim but different 
porosities V = 0.85, 0.74, 0.58, using three realizations per cluster (BA. 256. 1-3, BAM1.256. 1-3, and BAM2.256. 1-3) and 54 orientations per realization. 
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FIG. 7. — Applications to the debris disk around AU Mic. Scattering properties Sn(#)/Sn(0) and p(8) for 3 realizations of N = 256 BAM2 clusters 
(BAM2.256. 1-3), 54 orientations per realization, with monomer sizes oq ranging from 200A to 800A (R from 0.169/mi to 0.678/^m). Dashed lines show 
scattering properties obtained by averaging over a size distribution dn/ dR oc ff -3 5 running from R mm = 0.127/xm to i?max = 0.551 fim (see text). Solid curves 
show phase function and polarization inferred by Graham et al. 12007) for an assumed angular dependence given by eq. d4!5t . The adopted size distribution 
provides a good fit using clusters with porosity V 0.6 and sizes R 0.13—0.55 /im. For comparison, we also show results for clusters with oq = 700, 800 A, 
which exhibit the negative polarization branch observed in comets. Scattering properties available at http://www.astro.princeton.edu/~draine/SDJ09.html . 
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FIG. 8. — Scattering properties for aggregates resembling cometary dust (see text). For each composition we present results for N = 1024 BAM2 clusters (3 
realizations, BAM2. 1024. 1-3, 54 orientations per realization) with monomer sizes oq = 0.08 and 0.1 fim (R = 1.12/zm and 1.40/xm). Dashed lines are for a size 
distribution dn/dR <x R~ 3 - 5 for 0.98 fim < R < 1.54/xm. Scattering properties available at http://www.astro.princeton.edu/~draine/SDJ09.html . 



